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BeppoSAX Observations of the Power and Energy Spectral 
Evolution in the Black Hole Candidate XTE J 1650-500 

Enrico Montanari^, Lev Titarchuk^'^ and Filippo Frontera^ 

ABSTRACT 

We study the time variability and spectral evolution of the Black Hole Can- 
didate source XTE J1650-500 using the BeppoSAX wide energy range (0.12-200 
keV) observations performed during the 2001 X-ray outburst. The source evolves 
from a low/hard state (LHS) toward a high/soft state (HSS). In all states the 
emergent photon spectrum is described by the sum of Comptonization and soft 
(disk) blackbody components. In the LHS, the Comptonization component dom- 
inates in the resulting spectrum. On the other hand, during the HSS observed 
by BeppoSAX the soft (disk) component is already dominant. In this state the 
Comptonization part of the spectrum is much softer than that in the LHS (pho- 
ton index F is ~ 2.4 in the HSS vs. F ~ 1.7 in the LHS). In the BeppoSAX 
data we find a strong signature of the index saturation with the mass accretion 
rate which can be considered as an observational evidence of the converging flow 
(black hole) in XTE J1650-500. We derive power spectra (PS) of the source 
time variability in different spectral states as a function of energy band. When 
the source undergoes a transition to softer states, the PS as a whole is shifted 
to higher frequencies which can be interpreted as a contraction of the Comp- 
ton cloud during hard-soft spectral evolution. It is worthwhile to emphasize a 
detection of a strong low-frequency red noise component in the HSS PS which 
can be considered a signature of the presence of the strong extended disk in the 
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HSS. Also as a result of our data analysis, we find a very weak sign of Kq, line 
appearance in this BeppoSAX data set. This finding does not confirm previous 
claims by Miniutti et al. on the presence of a broad and strongly relativistic iron 
emission line in this particular set of BeppoSAX data. 

Subject headings: accretion, accretion disks — black hole physics — stars: individ- 
ual (XTE J1650-500): radiation mechanisms: nonthermal — physical data and 
processes 



Introduction 



XTE J1 650-500 was discovered by the Rossi X-Ray Timing Explorer {RXTE) on 2001 
September 5 (IRemillard et al. 11200 ll ) and subsequently reached a peak X-ray intensity of 0.5 
crab. Subsequent observations established that XTE J1650-500 is a stron g black hole candi- 



date based on its X-ray spectrum and variability in the X-ray light curve (IMarkwardt et al 



2001 



2003 



Rossi et al. 



Revnivtsev fc Sunvaev Il200ll : IWijnands et al. Il200ll : iHoman et al. 



IS; 



Kalemci et al. 



20041 ). Th e radio counte r part w as discovered with the Australia Telescope 



Compact Array (ATCA) by iGroot et al. I (120011 ). Furt her radio observatio ns sampled the 



behav ior of XTE J1650-500 along all its X-ray states (jCorbel et al. 1120041 ). lOrosz et al. 



(120041 ) found the period of the binary system to be 7.63 ±0.02 h, and the mass function [low 
limit of black hole (BH) mass] to be 2.7 ± 0.6 M©. 

We define the spectral state hardness according to the peak position in the EF(E)- 
diagram related to X-ray spectrum F(E) (in units of erg keV~^s~^) where E is the photon 



energy [see e.g. Figure 3 in lGrove et al. I (Il998l )]. Following this definition, we call a low/hard 
state (LHS) the spectral state in which the maximum of the EF{E)— diagram is located at 
tens of keV, while we call a high/soft state (HSS) the state where the EE {E)-maximnm 
is located at a few keV. Moreover, we call an intermediate state (IS) the state in which 
the source undergoes a spectral transition between the low/hard and high/soft states. We 
observed XTE J1650-500 three times with BeppoSAX during the 2001 outburst (September 
11-12, 21-23, October 3-4, see Table [1]), during which the source evolves from the LHS 
toward the HSS. The observed outburst is a perfect case to study a source transition from 
the LHS to the HSS. 

The long-term monitoring of Cyg X-1 has revealed not only two distinct energy spectra 
and transitions between them but it has also established that the corresponding time variabil- 



i ty po wer spectra (PS) evolve along with the photon spect r a [see IShaposhnikov fc Titarchuk 



(120061 ) ■ hereafter ST06, and IShaposhnikov fc Titarchuk I (120071 ). hereafter ST07]. Recently 
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the evolution of the power density spectra throughout the diffe rent X-ray states of BH X- 



ray bi naries has been also extensively studied in many BHs by iKlein-Wolt &: van der Klis 



fl2008[ ). hereafter KVK08. 



The persistence of the simultaneous PS and X-ray energy spectrum evolution suggests 
that the underlying physical process and conditions which give rise to the PS are tied to 
the corona; and, furthermore, that this process varies in a well defined manner as the source 
progresses from one spectral state to another. 

Moreover the same evolutions are seen in so many galactic X-ray binary BH sources 
(see, for example, detailed discussion of these evolutions in Cyg X-1, GRO 1655-40 and 
GRS 1915+ 105 in ST07) which vary widely in both luminosity (presumably with mass 
accretion rate) and state. This fact suggests that the physical conditions controlling the 
photon spectral and the PS evolutions are similar. In particular, quasiperiodic oscillation 
(QPO) frequencies are characteristics of the power spectra for these sources. Given that 
correlations between low and high frequency QPOs and between low f requency QPOs and 



break frequencies iri their PS have been found for a number of them (see lPsaltis et al. Ill999 



Belloni et al. II2002I ). one can suggest that QPO phenomenon can be an universal property 



of all accreting compact systems. 

In particular, the PS of black hole binaries in hard states is dominated by a component, 
which has a specific shape roughly described by a broken power- law (bknpl). The low- 
frequency part of the BKNPL is mostly flat, while the higher frequency part above the "break" 
frequency Uhr is a power law z/^" which index a varies between 1 and 2. In fact, the P S 



of BHs are more complicated than just a broken power law [see iTitarchuk et al. I (120071 ). 



Titarchuk fc Shaposhnikov I (120081 ). hereafter TSA07, TS08 respectively, and KVK08 for 



comprehensive study of the BH power density spectra in different X-ray spectral states]. 

It is also well established that the fractional root-mean-square (rms) variability in a 
source light curve decreases while QPO low-frequency (and Ubr) increases as the source 
evolves from the LH state to the HS state [see e.g. TS08, KVK08]. 



Lyubarskii I (119971 ) was the first to suggest a model for this time variability production 
in the accretion powered X-ray sources. He considered small amplitude local fluctuations 
in the accretion rate at each radius, caused by small amplitude variations in the viscosity, 
and then studied the effect of these fluctuations on the accretion rate at the inner disc edge. 
His linear calculations show that if the characteristic time-scale of the viscosity variations 
is everywhere comparable to the viscous (inflow) time-scale, and if the amplitude of the 
variations is independent of radius, then the power spectrum of luminosity fluctuations is 
a power-law of index a = 1, namely cx u^^. If the amplitude of the variations increases 
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with radius, the slope of the power spectrum of the luminosity variations a is steeper than 
1. Ly ubarskii pointed out that he h ad no physical model for the cause of such fluctua- 



tions. lUttley. McHardy &: Vaughan I (120051 ) pointed out that rms-flux relation is naturally 
explained in the framework developed by Lyubarskii. 

TSA07 formulated and solved the problem of local driving perturbation diffusion in a 
"disk-like" configuration (which can be either a geometric thin Keplerian accretion disk or a 
Compton cloud). The problem of the diffusive propagation of the spatially distributed high- 
frequency perturbations is formulated as a problem in terms of the diffusion equation for 
the surface density perturbations. This equation is combined with the appropriate bound- 
ary conditions. The formulation of this problem and its solution are general and classical. 
The emergent power spectrum as a result of the perturbation diffusion in a given bounded 
configuration (either Compton cloud or Keplerian disk) has a white-red noise (WRN) shape 
which can be approximated by a broken power-law (bknpl) with two indices and 
respectively. The parameters of the WRN PS are a timescale of the diffusion propagation of 
the local perturbations to and the power-law index of the viscosity distribution over radius. 

Note that the shapes of X-ray photon spectra of black hole (BH) and neutron star 
(NS) are generic: they consist of Comptonized blackbody-like (bb) components. These 
Comptonization components are seen as a power-law (pl) at energies higher than the BB 
characteristic energies. Sometimes one needs an exponential rollover in order to terminate 
the PL component at high energies. The main difference is that NS spectra are usually 
softer for the same state. Another difference in these spectra is that in the NS case there 
are two Comptonized BB components (not one as in BHC spectra) which can be related to 
the emission from the disk (BB color temperature about 1 keV) an d NS surface (BB color 



temperature about 2.5 keV [see iTitarchuk. fc Shaposhnikov I (120051 )]. 



Titarchuk et al. I (119971 ). hereafter TMK97, introduced the Generic Comptonization 
model (bmc model in XSPEC) that takes into account the dynamical Comptonization (con- 
verging inflow, expected in the vicinity of the central object) along with the thermal Comp- 
tonization. Note that this model can be applied to fit the observed spectra of NS and BHC 
sources. 

In this Paper we present the results of our evolution study of the energy and power 
spectra in XTE J1650-500 using the BeppoSAX observations and using the aforementioned 
theoretical considerations of photon and power spectral formations. We find that this evo- 
lution is similar to that observed from Cyg X-1 (see e.g. ST06 and TSA07). 



Miniutti et al. I (120041 ). hereafter M04, reported spectral results from three BeppoSAX 



observations of XTE J1650-500 during its 2001/2002 outburst. By performing the analysis 
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in the energy range from 1.5 to 60 keV for the TOO-1 spectrum and in the range from 1.5 
to 200 keV for the other two spectra, they reported the presence of a broad and strongly 
relativistic Fe hne with an EW of about 200 eV, depending on TOO. We reanalyze the 
same BeppoSAX data set adopting the entire BeppoSAX energy range (0.12—200 keV) in 
order to check the presence of a strong K„ line in the data. We find a very weak sign of 
line appearance. 

Results of our detailed spectral analysis are the subject of our next paper (Montanari 
et al. 2008 in preparation). Here we present results on the simultaneous evolution of the 
spectral and temporal properties and compare them with that in other black hole (BH) 
sources (see ST06 and ST07). 

Our approach is two-fold: On one side we tried to fit all photon and power spectra with 
the same generic simple model, in order to compare them directly. On the other hand, we 
study the source in the framework of a physical model in the attempt to find a self-consistent 
scenario that can help us to understand the Physics of the spectral and timing evolution that 
we observe. 

In §2 we briefly describe the BeppoSAX data from XTE J1650-500. Details of our 
spectral fitting are presented in §3. We discuss an observational evidence of a converging 
flow in §3.1. We present the results of our analysis of Kq, appearance in the data in §3.2. 
Description of XTE J1650-500 power spectra and QPO identification are presented in §4 
where we also discuss a correlation of PS vs spectral state (photon index) revealed in XTE 
J1650-500. Discussion and conclusions follow in §5. 



Observations and data reduction 



XTE J1650-500 was observed on September 11-12, 21-23 and October 3-4 2001 during 
three Target of Oppor tunity Observations (TOOs) by the Narrow Field Instruments (NFIs) 
on board BeppoSAX (IBoella et al. Ill997al ). Table [T] reports the observation log along with 
the on-source expo sure times. The NFI s include a Low-Energy Concentrator Spectrometer 



(LEGS, 0.1-10 keV: IParmar et al. 



(MEGS, 1.5-10 keV; iBoella et al. 
Gounter (HPGSPG , 4-120 keV; 



19971). three Medium-Energy Goncentrator Spectrometer 



1997br). a H i gh-P ressure Gas Scintillation Proportional 



Manzo et al. 



19971 ). and a Phoswich Detection System 



(PDS, 15-200 keV; iFrontera et al. 1119971 ). The SAXDAS data analysis package is used for 
processing data. 



For each instrument we perform the spectral and temporal analysis in the energy range 
in which the response function is well determined; given the high statistics of the source the 
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energy range is 0.12-4.0 keV for the LEGS, 2.5-10 keV for the MEGS, 8-30 keV for the 
HPGSPG, and 15-200 keV for PDS. Note that, at the given high count rate, the available 
LEGS response function is not reliable, therefore the SAX/LEGS Matrix Generation Program 
"lemat" is used to obtain the appropriate response function of the instrument. MEGS unit 
1 was not operative. Spectra from MEGS unit 2 and 3 are merged. The spectra so obtained 
are rebinned taking into account the energy resolution of the instruments in order to get 
independent data points. 

Previously M04 analyzed the same set of BeppoSAX data integrating the source spectra 
over the entire duration of each of the three observations. However, this integration time 
could provide distorted spectra, given that the source shows a significant spectral variability 
with time on shorter time scales. This effect was first revealed studying the time variability 
of three different hardness ratios i?[(4 - 10)keV/(2 - 10)keV], R[{10 - 30)keV/(2 - 10)keV)], 
R[{20 - 60)keV/(2 - 10)keV] with bins of 200 s integration time. The results, that will be 
reported in detail elsewhere (Montanari et al. 2008, in preparation), show that, depending 
on the epoch, significant (up to 30%) spectral variability is present in time intervals ranging 
from 1 to 10 ks. Therefore, unlike M04, we integrated the spectra over the longest time 
intervals compatible with a negligible (less than 5%) source spectral variability. 

Namely, in the TOO-1 the spectra were divided in 12 time intervals and numerated from 
1101 to 1112, in the TOO-2 they were divided in 19 intervals (from 1201 to 1219), while in 
the TOO-3 they were divided in 8 intervals (from 1301 to 1308). Note that LEGS data are 
not available for five intervals of TOO-2. 

In Figured] we display the source light curve related to the period from September 11 to 
October 4, 2001 for the 2-3 keV and 20-60 keV energy bands. The integration time of each 
bin is that of the IDs intervals used for spectral analysis. In Figure [T] one can clearly see 
an overall anti-correlation between the two light curves related to these two energy bands 
during TOOl and TOO-2. Note that the transition to a softer spectrum becomes prominent 
in TOO-2 and continues during TOO-3 observations. For TOO-2 (see Fig. [1]), a variable 
rate is evident at low energies, unlike an almost stable count rate at high energies. This 
different time behaviour at different energies can be also seen from Fig. [2] where we compare 
the energy behavior of MEGS mean count rates in the time intervals 1202 and 1207. 



We model the energy spectra using XSPEG 11.2.0 software package (lArnaud 1 119961 ). 
A systematic error of 1% is added in order to take into account unavoidable uncertainties 
in the response functions. In the multi-instrument fits we leave free to vary the relative 
normalization of LEGS, HPGSPG and PDS with respect to MEGS and the derived unfolded 
spectra are renormalized to the MEGS level for a clarity of display. 
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For the temporal analysis we use standard fast Fourier techniques. Source PSs are 
normalized to rms units using subtraction of the Poissonian noise level. Dead time corrections 
have been taken into account according to Ivan der Klis I (119891 ). The binning time for the 
analysis is a few ms. PSs are calculated using data stretches of ~ 100 s and then averaged. In 
this way we obtain the reliable PSs in a frequency range from 0.01 to 100 Hz. Unfortunately 
we cannot provide much longer time stretches to extend the frequency pass band to lower 
values due to frequent interruptions in the data rate. We present only data points which are 
statistically significant (see more details in §4). 

The uncertainties of the model parameters are given at 90% confidence level, while those 
shown in all figures are la errors. 



3. X-ray energy continuum spectra 

In our analysis, all derived spectra (39) cover the 0.12-200 keV energy band, except 5 
intervals in which the LEGS data are not available. In this case the energy band was from 
2.5 to 200 keV. 

To describe t he continuum spectruni we use the Bulk Motion Comptonization model 



(bmc in XSPEC) (iTitarchuk et al. 1119971 ) which is a generic Comptonization model. This 
model can be used if the photon energy is less than the mean electron energy of the Compton 
cloud Eav The choice of this particular theoretical model is provided by the robust nature 
of the BMC model for different spectral states and independence of the specific type of 
Comptonization scenario involved. The BMC model spectrum is the sum of a blackbody 
(BB) component (which is the disk radiation directly seen by the observer) and the fraction 
of a BB component Comptonized in the corona. 

The model has four parameters: the color temperature kT^f, of a thermal photon spec- 
trum (blackbody, BB), the energy spectral index a = F — 1, where F is the PL photon 
index, the parameter A, which is related to the weight [A/{1 + A)] of the Comptonization 
component, and the normalization of the BB component Cn- The BMC model is valid for 
the general case of Comptonization when both bulk and thermal motion are included. 

For the thermal Comptonization, Eav is related to electron temperature Eav ~ 2kTf,. 
When the bulk motion Comptonization is dominant, Eav is related to the bulk kinetic en- 
ergy of the electrons EavtmeC^ ■ In the thermal Comptonization case, for energies less than 
Eav, COMPTT0 and BMC models are very similar. The thermal Comptonization and the 



^COMPTT model was introduced by IXitarchuk I (|l994l ). IXitarchuk fc Lvubarskii I (|l995( ) and 
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dynamical (bulk motion) Comptonization are presumably responsible for the spectral for- 
mation in the hard state and the soft state, respectively. Therefore, one needs a general 
model such as the BMC that describes the spectral shape regardless of the specific type of 
Comptonization. Although the thermal Comptonization model can properly fit the spectral 
shape of the observed spectra for all spectral states, it would give ph ysically unreasonable 



values of the best-fit parameters for the soft state spectra. Notably, IWilms et al. I (120061 ) 
found optical depth r of the Compton Cloud inferred from the COMPTT model that signif- 
icantly decreases towards the HSS. It is very difficult (in the framework of any reasonable 
physical model) to justify this tendency of r to decrease when the mass accretion rate in- 
creases during the hard-to-soft state transition (see e.g. ST06). However the spectral index 
a inferred using the best-fit COMPTT parameters r and fcTe is a physical characteristic of 
the Comptonization process. In fact, the spectra l index a is the reciprocal of the Comp- 



tonization parameter Y [see e.g. iBradshaw et al. I (120071 )] and it is independent of any type 



and model of Comptonization (thermal or bulk). 

It is important to emphasize the principal difference and similarity in using blackbody- 
like spectra (bbody in XSPEC) and mu l ticolor disk spectra (diskbb in XSPEC) as soft 



photon spectra. In fact, iBorozdin et al. I (Il999l ). hereafter BRT99, demonstrated that in a 



limited energy range (for example RXTE/VCA) the diskbb is fit by a BB spectrum. The 
color BB temperature is the best-fit parameter which is related to the effective radius of the 
emission area of the disk. It is evident that the best-fit values of the color temperature and 
effective radius depend on the energy range of the instrument (see BRT99). Any disk model 
requires the integral of the local spectra over radius. These local emergent spectra are not 
necessarily blackbody-like spectra and they are indeed formed as a result of the radiative 
transfer effects which depend on the assumption of the local density and temperature struc- 
ture. Moreover the specific diskBB model assumes a temperature distribution with the disk 
radius that is likely not the real distribution. Thus, taking all that into account, we have 
preferred to describe the thermal (BB-like) component of the source spectrum by adding 
another BB to that already included in the BMC model. 

Following XSPEC notation, the model used to describe our spectra is given by F{E) = 
WABS X {BBODY + BMC) x SMEDGE x HIGHECUT. Here WABS is the XSPEC 
photoelectic absorption model which takes into account a possible intrinsic photoelectic ab- 
sorption in the source and/or external absorption of X-ray source emission in the way towards 
the Earth. The BMC Comptonization model describes both direct (from an innermost part 
of the disk) and Comptonized photons (from the corona), the added bbody component 



Hua fc Titarchuk I (|l995h to describe the spectrum of thermal Comptonization in the whole energy range. 
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describes the direct soft photons belonging to the disk part at a greater distance from the 
central object which spectrum can be only slightly affected by upscattering in the corona. 
The SMEDGE (smeared K-edge) model represents a possible in teraction of the photons orig - 



inated in the inner part of the X-ray source with a wind (see iLaurent fc Titarchuk 1120071 ). 
The HIGHECUT (high-energy cutoff) model enables to take into account the curvature effect 
of the Comptonized photon spectrum at high energies (that is above 100 keV in the LHS 
spectra) due to the recoil effect which is neglected by BMC. highecut is found to be needed 
to describe the LHS energy spectrum, while, for the HSS, no high energy cutoff is needed. 

Most of the 39 spectra are found to be very well fit by the assumed input model. Only 
for 9 spectra the model has a null hypothesis probability less than 0.05, with 3 of them with 
null hypothesis probability less than 0.01, the smallest value being 3.7 x 10~^. 

The best-fit parameters of 6 typical spectra extracted from the entire set are shown in 
Table [2l As can be seen, the temperature kTstb of the soft blackbody varies from 0.18 to 
0.36 keV, while that of the BMC, kT^b, varies in the range from 0.39 to 0.56 keV, depending 
on the spectral state. For the same time intervals of Table [2], we show in Fig. [3] the best 
fit EF{E) spectra in keV cm~^ s~^ units. In the LHS spectra 1105 and 1111 the relevant 
Comptonization contribution is related to the thermal Comptonization of soft (disk) photons 
off hot electrons of Compton cloud (corona). Instead, in the HSS spectra (see spectra 1215, 
1216, 1303, and 1306 in Fig. [3]), the blackbody emission (both EEs), presumably related to 
the intrinsic disk spectrum, is dominant in the resulting spectrum, while the hard (steep pl) 
photons, probably formed in the relatively cold (a few keV temperature) compact region, 
are only a small fraction of the total flux. 

As we noted above, the use of the thermal Comptonization model (COMPTT) is phys- 
ically justified for the fitting of the LHS spectra only. For comparison with the results 
reported in Table [21 for time intervals 1105 and 1111 we implemented a fit using a model 
WABS X {BB + COMPTT) x SMEDGE (see Table E]). In fact, the fit quality obtained 
using the COMPTT model are not so good as that in the case of BMC [compare the x^/ (dof ) 
values in Table [2] and [3] presumably because a modest converging flow contribution to the 
overall Compton upscattering is already present in TOOL In general, the worse quality fits 
introduce additional systematic errors in the best-fit model parameters. In particular, the 
photon index, obtained using COMPTT best-fit optical depth r and electron temperature 



kTp, an d using the formula for the index related to r and kTg [see iTitarchuk fc Lyubarskij 



( 1l995l )]. is 0.4 higher than that inferred using BMC. 



In the sections 4-5 we present the arguments for the compactness of the steep PL emission 
region which is presumably a converging flow streaming towards the innermost part of the 
EH. 
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3.1. Observational evidence of the converging flow (black hole) in XTE 

J1650-500 



It is important to emphasize once again (see above) th at the spectral index a 



=r-ii 



IS 



the reciprocal of the Comptonization parameter Y [see e.g. iBradshaw et al. I (120071 )] and it is 
independent of the Comptonization model (thermal or bulk). The parameter Y = Ngj] is the 
mean number of scatterings Ngc times the average fractional energy change per scattering rj. 
The number of scatterings of Comptonized (upscattered) photons in the converging flow (CF) 
Nsc is proportional to the optical depth r, given that in the converging flow the photons can 
be only eflectively upscattered in the direction of the bulk motion, and thus Ng^. oc R/l = r, 
where R is the CF characteristic radius and / is the photon free path. It is also worth 
noting that, in the converging flow, r is also equal to the d imensionless mass accre tion 



rate rh = M/M^d where M^d is deflned as M^d = L^d/c^ [see iTitarchuk et al. I (119971 ) for 



more details]. On the other hand, the change of the CF fractional en ergy per scattering rj is 



i nvers ely proportional r (r^ oc 1/r) when r ^ 1 [see for more details in lLaurent fc Titarchuk 
(|2007h ]. 

Consequently Y = Ngcf] and thus F = 1/Y + 1 saturate to a constant value when r 
(or rh) increases. Observationally this index saturation can be revealed if one can flnd a 
correlation between F and BMC normalization Cn = L^g/dl^ where L 



39 



Lrf/(10^^ ergs" 



is the luminosity of the seed photon source (disk) in units of 10 erg s~ and dio is the 
distance to source in units of 10 kpc. This correlation of the index F with = L^g/dff^ 
implies the correlation of F vs m given that as a disk luminosity is proportional to m, 
namely Ld oc Seffmrh, where m = M/Mq is the black hole mass in units of solar mass 
and Sejf is the efliciency o f the gravitational energy release in the accretion disk [see e.g. 
Shakura fc Sunyaev I (119731 )]. In principle, the efliciency Seff can depend on the inner radius 
of the disk Rin^, i.e. indirectly depen ds on m. However, in HSS, R in d presumably reaches 
its lowest limit of about 3i?s [see e.g. IShrader fc Titarchuk I (119991 )] and thus, in HSS, the 
index vs. Cn correlation is equivalent to the correlation between index and m. This implies 
that in HSS the index saturation with Cat is equivalent to the index saturation with fn. 

In Figure m we present an observational evidence of the index saturation with L^g/dl^, 
and thus with rh, which can be considered as an observational signature of a converging flow 
(black hole) m XTE J 1650-500. 



3.2. On iron fluorescence line determination 



Limiting our considerations to the BeppoSAX TOOs we confirm tlie broad excess at 
~ 6.4 keV found by M04 in all of the three observations if we adopt, as M04 did, the XSPEC 
reflection model PEXRiv. However, the excess is present only in the first TOO if we use 
a two-phase accretion disk corona model which also includes the reflection effect (XSPEC 
COMPPS Model) (Montanari et al. 2008, in preparation). It is worth noting that when we 
use COMPPS model we find a strange effect (unexpected in the scenario described by this 
model) that the iron line (thought to be due to reflection from the disk) is present when 
the reflection parameter is low (in TOOl) and it is c ompletely absent when the reflection 



parameter becomes larger (in T002 and TOOS) (see iMontanari et al. II2004I ). Also as we 
have already mentioned above, in the framework of the COMPPS model there is a puzzling 
effect of significant decrease of the optical depth of the Compton cloud when the source 



evolves to the HSS (IMontanari et al. II2004I ). 



On the other hand, using the input continuum model described in Section [3] and the 
entire BeppoSAX spectral energy range we find a very marginal evidence for the presence of 
the Fe fluorescence line in the data. We illustrate our findings in Figure [51 

Adopting the 0.12-200 keV energy range, the fit of the continuum model WABS x 
{BB + BMC) X SMEDGE x HIGHENCUT to the spectrum 1108, used as an example, 
is shown in the panel 2 (starting from the top). As can be seen, the above continuum model 
gives a very satisfactory description of the data (x^/dof = 150/164), with no need of adding 
a Fe Kq, line model. Instead the SMEDGE model appears crucial for a good fit. The fit 
of the above continuum model with no SMEDGE to the same spectrum (see top panel) is 
highly unacceptable {x^/dof = 430/165). The apparent sinusoidal pattern of the residuals 
to the model is a clear sign of the bad fit quahty. 

Note the addition of a Kq, line (laor model) to the continuum model with no SMEDGE 
gives an unacceptable fit (x^/dof = 300/163) (see panel 3), but it becomes acceptable 
(x^/dof = 136/109, see panel 5), if we limit the spectrum energy bandwidth to 1.5-60 keV 
(the band used by M04) . However the model which is valid in 1.5-60 keV, is not valid in 
the entire BeppoSAX energy band. Namely the model WABS x {BB + BMC + LAOR) x 
HIGHENCUT, with the best-fit parameter values found for the 1.5-60 keV energy band, 
does not fit the data related to the entire BeppoSAX energy range from 0.12 to 200 keV 
(xVdof = 1315/164, see panel 6). 



^Note that a fit of the model WABS x {BB + BMC) x HIGHENCUT with no Laor component does 
not fit the data in this limited energy band (see panel 4). 
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We have also attempted to understand why M04, using the SMEDGE model as we 
do, obtained different results. Limiting the analysis to the spectrum 1108, assumed as an 
example, we actually find that, freezing the energy of the Edge to the best-fit value re- 
ported by M04 \ Epd,ge, = 8.3 keV b eing c onsistent with a very high ionization state of iron 



ion XXII), see iKallman et al. I (120041 )] gives the presence of a significant emission line. 



The line disappears when Eedge is left free to vary, finding, as a best-fit value, 7.2 keV. 

To conclude, iron line appearance strongly depends on the energy band used to 
describe the continuum emission. With our best-fit continuum model, that gives the best 
description of the source spectrum in the entire 0.1-200 keV BeppoSAX energy band, the 
line is not required. 

Although no significant emission line is evident in the data, adding the XSPEC LAOR 
model to our best-fit model (see Section 3) slightly improves the fits (Ax < 10 with a 
decrease of 4 dof for about 160 initial dof) for some observation intervals belonging to TOOl 
(LHS). The resulting line, of equivalent width ranging from few tens eV to about 100 eV, is 
almost symmetric with its center being situated in the energy range between 5 and 5.5 keV. 
For all of other intervals the normalization of the line goes to zero. 

Concerning the SMEDGE model used in the best fit description of our spectra, it is 
worth noting there is a relation between Thompson optical depth tt of the iron emission 



area and the SMEDGE maximum absorption optical depth Tmax at threshold [see iBasko 



([19781) and ST06]: 



/7.8keVV 

Tmax = {Y + Yo) ( " ^ ) (1) 

where Y and Yq (in units of the cosmic values) are the abundances of the elements with 
a charge Z < 26 and the iron abundance respectively, while Eth is the K-shell ionization 
threshold energy. In our case F ~ 1, Fq — 1, Eth = 7.1 keV, so 

Tmax - 2.65rT 

For instance, Tmax = 2, 3, 4, 5 correspond to tt = 0.75, 1.1, 1.5, 1.9 respectively. 

We find that, in the LHS, the Thomson optical depth tt is relatively low, i.e. about 
0.7 whereas, in the HSS, tt is about 2 (see Table [2] for the values of Tmax and Eq. [T]for the 
Tmax — Tt relation). Thus, while in the LHS the possible detection of the line is consistent 
with the value of the Thompson optical depth tt ~ 0.7 (see Table [2]), in the HSS, its non 
detection is consistent with the smearing due to a significant Thompson optical depth tt in 
the range between ~ 1.3 and ~ 2 (see above). 

One can ask a fair question on the relation between the K-edge and appearances in 
the data. In many models Ko-line formation is indeed strongly related to edge formation 
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[ see e .g. 
(53), 



Basko et al. I fll974r) ISasko I fllQTSh . iKallman et al. I fl2004f ). iLaming fc Titarchuk 



Laurent &: Titarchuk 



(120071 )]. However it is not always true that the detection of 
the K-edge in the data imphes the Kq,— hne detection. The physics of the K-edge and 
formation is the following: photons of energies higher than 7.1 keV ionize the iron K-shell 
and this ionization effect leads to emission of K^— photon (iron K^— fluorescence) with a 
probability of ljk (where ujk is K-yield), whereas with probability (1 — uk) the 
energy goes to on ionized electron. So the edge can be always there but the K^— photon is 
emitted with a probability less than 0.3 and in addition this probability is affected by the 
electron density of the surrounding medium (see more details in the aforementioned references 
on this subject). Even if the K^— photon is emitted it does not mean that this photon can 
be detected and seen by the Earth observer. These photons can be absorbed when they 
propagate through the medium (for example in the wind) and they can be also scattered off 
electrons and significantly change their energy, i.e. washed out from the resulting spectrum. 
A H these effects of line appeara nce and disapperance have been reproduced in simulations 
by iLaming fc Titarchuk I (120041 ) and iLaurent &: Titarchuk I (120071 ) . 



4. Power spectrum and timing-spectral correlation 

We succeeded to derive accurate power spectra (PS) of the source time variability thanks 
to the high source brightness during outburst. In particular, during TOO-1, when the source 
was mainly in the LHS, it was possible to obtain the PS in several energy intervals within 
the broad 0.1-200 keV BeppoSAX operational band. 



4.1. Phenomenological (Lorentzian) Model 

First we made the analysis of the power spectra using, as an input model, the superpo- 
sition of Lorentzian functions (Lorentzian model). We compared their characteristics with 
those identified by KVK08 for a number of BHs. A Lorentzian function is given by 

Liu) = K . (2) 

[v- + (r/2)2 

where z/q is the centroid frequency and P is the full width at half maximum of the function 
(FWHM). 



We used broad Lorentzians to fit the continuum spectrum and narrow Lorentzians to 
fit QPO features. In Figure E] {bottom panel) we illustrated the fitting results and the 



-14- 



characteristics of the PS components in terms of frequency x power diagram. Note the power 
spectrum unit is (rms/mean)^/Hz and thus the frequency xpower unit is (rms/mean)^. 
Below we presented a comparison of these characteristics with those found by KVK08. 

The continuum of all derived LHS PSs is found to be well described by the sum of 
two broad Lorentzians Li and L2. Following the KVK08 notation, we formally call them 
Lh and Lh, respectively. We report the best fit characteristic frequency of these broad 
Lorentzians in Table |H This frequency corresponds to the maximum of the vL{v) diagram 
^max = [i^o + (r/2)^]°'^. The presented PSs are results of integration over the entire duration 
of each TOO. 

Using Fig. [6] and Table H] one can see that, in addition to the two broad Lorentzians 
Lfc and L^, the 2-10 keV power-frequency diagram of the LHS (TOOl) requires two narrow 
Lorentzians: a low frequency (LF) QPO and a higher frequency QPO at vqpoi ~ 1-5 Hz 
and i'Qpo2 ~ 3 Hz, respectively. These two QPOs were introduced by KVK08 as LF and 
LF"*", respectively. Note that i'qpo2 ~ 3 Hz (or LF+) was found to be consistent with the 
second harmonics of i^qpoi (or LF). 

However, in the 0.1-2 keV energy channel, no evidence of vqpo2 is found and the 
characteristic frequency of the broad Lorentzian L/^, that describes the hump to the right of 
the QPO frequency (see blue line in the right bottom panel of Fig. [6]), is at 23 Hz. Note the 
L/i— hump is located at ~ 2 Hz in the 2-10 keV power-frequency diagram. 

In Table [5] we show the energy dependence of the LHS PS components in the broad en- 
ergy bands from 0.1-2 keV to 15-200 keV. The characteristic frequencies of L\, are consistent 
with each other up to 15 keV and the Lh characteristic frequencies are consistent with each 
other in the 2-4, 4-10 and 8-15 keV energy intervals. Us ing the Lh — Lh relation found (f or 



the energy band 2-25 keV) by Wijnands and van der Klis |Wijnands fc van der Klis I (119991 )]. 
updated in KVK08, we find that our medium energy data points are consistent with that 
relation. Also the PS narrow components Lgpoi and \^qpo2 clearly detected in the 4-10 keV 
and 8-15 keV bands satisfy the Wijnands and van der Klis L^p — Lb relation (see Fig. 10 in 
KVK08). 

It is worth noting that, in the 0.1-2 keV and 15-200 keV energy intervals, the Lh 
frequencies are an order of magnitude higher than those at medium energies. In fact, it is 
the first time that the high values of I'max (~ 23 Hz, > 50 Hz in the 0.1-2.0 keV and 15-200 
keV energy bands respectively), are discovered from this source. Their identification with 
the Lh introduced by, e.g., by KVK08, is not obvious. Unfortunately, there are no data in the 
literature to statistically study a possible correlation between Lb and Lh at these energies, 
unlike what can be done at medium energies with RXTE data (see, e.g., KVK08). The Lh 
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components in the 0.1-2 keV and 15-200 keV energy bands can be related to PS physical 
components which are different from that in the medium energy bands (see details of the 
physical interpretation in section 4.2). 

In the HSS, the PS continuum is still found to be well described by the sum of two broad 
Lorentzians, but in this case (see Table H] and Fig. [6]), L/^ is no more visible, while a very 
low frequency noise (VLFN) emerges. It is important to emphasize that the characteristic 
frequency of L;, increases from the LHS to the HSS. Our finding is in agreement with the fact 
that the Lf, characteristic frequency is shifted to higher values during the source evolution 
from the LHS to the HSS [see also Fig. 4 in KVK08]. We also confirm the KVK08 result 
that LFVN component occurs in the HSS, but is absent in LHS [compare our Fig. ^{bottom 
panels) with panels a and e in Fig. 4 of KVK08]. 

Furthermore we find that the Q- value, defined as Q = I'o/T, is in the range of 0-0.5 
for Lb, Lh components and that is in the range of 4-6 for PS narrow components Lqpq- As 
one can see from Fig. [6] that rms/mean values are 10% and 30% for Lb and L^ of the LHS 
power spectra components respectively in 0.1-2 keV. The corresponding rms/mean values 
for the LHS Lf, and in 2-10 keV band are about 10% (see left bottom panel of Fig. [H]). 
However the variability power of the PS components drastically decreases in the HSS, and 
their rms- values are less than 1%. It is worth noting that these ranges of Q and rms values 
are consistent with that found by KVK08 (see Table 1 there). 

4.2. Physical (diffusion) model 

The analysis of the PSs was performed using a simplified version of the diffusion model 
(see TSA07 and TS08) in which the PS continuum shape at frequencies below the driving 
frequency can be approximated by a broken power-law (bknpl). The results of the model 
fitting to the data are shown in Tables [MHl We also present the model fitting results in 
Figs. EE 

In particular, in Fig. [7] and Table [6] we show the energy dependence of the LHS PSs and 
the variability power, that decreases with energy (mainly above 10 keV), while it preserves 
the self-similar shape for energies higher than 2 keV. These results obtained for the LHS 
have relevant implications. The self-similar shape of the > 2 keV spectra implies that 
photons of these energies are produced (upscattered) in the same geometrical configuration. 
Following TSA07, the break frequency found in the LHS PSs is related to a diffusion time of 
perturbation propagation while the QPO low frequency is an eigenfrequency of the volume 
(magnetoacoustic) oscillation of the medium (in our case it is a Compton cloud). Our results 
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show that the values of these PS characteristics are similar in all PSs for photon energies 
higher than 2 keV. 

As we have already demonstrated in Section 3, the energy spectral behaviour (Fig. [3]), 
even within TOO-1, shows a slight softening from the time interval 1105 to 1111. In order 
to see the effect of this softening in the PS, we have compared the 4-10 keV PS derived in 
the time interval from 1101 to 1107 with that related to the time interval from 1109 to 1112 
in the 4-10 keV energy band. Note that the 4-10 keV band was chosen because we find the 
best signal to noise ratio in this interval. 

The result is presented in Table [7] and in Fig. [81 The main effect shown in Fig. [8] is the 
shift to higher frequencies of the PS which correlates with the photon spectral softening (see 
Fig. [3]). Independently of the continuum adopted model the LF QPO centroid frequency 
increases while the continuum noise remains stable in its shape. For the diffusion model 
continuum, as can be seen from Table[7l the LF QPO frequency increases from 1.47±0.03 Hz 
to 1.66 ± 0.08 Hz (uncertainties at 90% confidence level). 

In order to illustrate the evolution of the power spectra from TOO-1 to TOO-2, when 
the energy spectra undergo their transition from the LHS to the HSS (see 1215 and 1216 
spectra in Fig. [3]), we compared the LHS PS with the HSS PS in two energy bands where 
the signal is highest: 2-10 keV and 0.1-2 keV (see top left and top right panels of Fig. [6l 
respectively). As can be seen from this Figure, the relative power of the HSS PS (with 
respect to the LHS PS) and its shape strongly depend on the energy band. 

We find (see Table [8]) that the 0.1-2 keV HSS PS is well described by the sum of a 
PL plus a BKNPL, with no evidence of QPOs in contrast to the 0.1-2 keV LHS PS which is 
described (see Table E]) by a bknpl plus a high (z/q ~ 14 Hz) frequency broad Lorentzian 
(FWHM~ 40 Hz). On the other hand the 2-10 keV HSS PS (see Table El and Fig. E> can 
be fit with a PL with plus a constant, while the 2-10 keV LHS PS is fit with a bknpl plus 
two narrow Lorentzians, one of which is consistent to be the second harmonics of the other 
at z/q ~ 1.5 Hz. 

These results confirm what was previously found by iHoman et al. ] J2003[l using RXTE 



data, that t here are prominent QPO features in XTE J1650-500 PSs [see also 



mm and lWiinands et al. I fl200lh ]. 



Revnivtsev fc Sunyaev 



In the LHS, the Comptonization component is dominant and thus we should detect 
the X-ray emission along with its variability related to the Comptonizing region (Compton 
cloud) only. Indeed in all of the LHS PSs shown in Fig. [6l the VLFN component is not seen. 
Its absence is in agreement with the low intensity level of the direct disk blackbody component 
in the LHS energy spectra (see 1105 and 1111 spectra in Fig. [3|). It is important to emphasize 
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that, in the LHS, the typical values of the characteristic frequencies of the VLFN and bknpl 
components are of the order of 10"'' — 10~^ Hz and 0.1 — 1 Hz, respectively (see TSA07 for 
Cyg X-1 case and KVK08). Thus, in this spectral state, the disk variability can be clearly 
detected only if frequencies lower than 10~^ Hz corresponding to photon energies less than 
~ 0.1 keV (see TSA07). However this is not the case for BeppoSAX . 

On the other hand in the HSS PSs the VLFN component appears (with the power 
depending on the energy range) along with the bknpl component. Note that this latter 
component is barely visible in the 2-10 keV PS (see Fig. E]). This result is in agreement 
with a low contribution of the high energy tail in the HSS energy spectra (see Fig. [3]). 
Following TSA07, we can interpret the VLFN component observed in the HSS as related to 
the timing response of the extended disk, while bknpl or (in terms of Lorentzian model) 
is associated with mass accretion rate perturbations in the Compton cloud (see the HSS PS 
in the right panels of Fig. [6]) . 

We find that in the 0.1-2 keV band not only Vb,max (using the Lorentzian model) but 
also V}yr (using the diffusion model) increase from TOO-1 (LHS) to TOO-2 (HSS) by a factor 
60 (z/f, max increases from 0.591^q'53 to 36 ± 3 Hz, while p^jr increases from 0.28 ± 0.03 Hz 
to ~ 17^2 Hz). This high shift of the related frequencies is in agreement with the corona 
(Compton cloud) contraction during the LHS-HSS transition (see TSA07). Searching in the 
literature, we found that in Cyg X-1 PSs the break frequency changes from ~ 0.03 Hz to 
~ 18 Hz value when photon index varies from 1 .5 to 2.1 (see Fig. 8 in ST06). Also we 
should point out that IWijnands fc van der Klis I Jl999h ( see Fig. 2) reported the range for 
the break frequency in BHs is 0.03 — 30 Hz. 

The ~ 14 Hz hump clearly seen in the 0.1-2 keV LHS PS (see Fig. [6] and Tables HI [6]) 
disappears at higher energies presumably because of the lack of statistics (compare with Fig. 
[7]). TSA07 and TS08 give an explanation of this hump in the high frequency tail of the power 
spectrum of Cygnus X-1, in terms of driving frequency of the perturbation oscillations of 
the Compton cloud. 



5. Discussion and Conclusions 

In this Paper we present the results of simultaneous spectral and timing analysis of the 
BHC XTE J1650-500 observed by BeppoSAX. 

The energy spectra are well described by the sum of a Comptonization and soft black- 
body components along with a smeared edge. In particular, we argue that the iron line 
appearance strongly depends on the energy band used to describe the continuum emission. 
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Using both the high energy (>20 keV) and the low energy (<1.5 keV) bands, which are 
both crucial to describe the real continuum level, the presence of the iron line feature is 
not mandatorily required. Thus we conclude that the line feature is not so prominent in 
BeppoSAX data as initially reported by M04. One may also attribute the difference of our 
conclusion regarding the presence of the Kq,— iron line in the BeppoSAX data and that by 
M04 to a different treatment of systematic errors and different procedure of data integration 
over time. 

We find the clear evidence of an evolution of the spectra from the LHS to the HSS (see 
Fig. [3]). In the LHS the Comptonization component dominates, in contrast to the HSS when 
the blackbody component dominates the emergent spectrum. Moreover in the BeppoSAX 
data (see Fig. H]) we find a strong signature of the index saturation with BMC normalization 
L^g/dlf), which is, in fact, proportional to the disk mass accretion rate. This index saturation 
vs. mass accretion rate can be co nsidered as an observa t ional evide nce of the converging 



flow ( black hole) in XTE J1650-500 Titarchuk fc Zannias I (119981 ) and iTitarchuk. fc Fiorito 



(120041 ) hereafer TF04]. 



The relative contribution of the Comptonization component (in the resulting spectrum) 
decreases when the source undergoes the hard-to-soft spectral transition (see Fig. [3]). We 
interpret this as a result of contraction of the Compt onization region (Compton cloud) during 



this spectral evolution [see ITitarchuk et al. I (119981 ). hereafter TLM98, and TF04] 



We can infer this contraction effect using our best-fit model parameters. One of the 
parameters of the BMC model is the normalization of the blackbody (BB) component Cn 
which is directly related to the effective area All^ of the BB emission [All^ oc CnHT^i^ (1 + 
A))]. If we multiply this effective area with the fraction of Comptonized component / = 
A/{l+A), we obtain the effective area of the high energy X-ray emission region (the Compton 
cloud). Using the best-fit values for the LHS and HSS spectra we find that the area of the 
Compton cloud decreases by a factor five during the transition from the LHS to the HSS 
(see Fig. W)- 

Furthermore in Fig. [TU] we demonstrate how the Comptonization fraction / decreases 
when the photon index increases. This correlation presents one more argument for a Comp- 
ton cloud contraction during LHS-HSS evolution. While in §4 we argued this contraction 
effect using the power spectrum analysis only, here we find more arguments for this contrac- 
tion using the energy spectrum analysis. 

Actually the power spectra follow the energy spectral evolution. The PSs related to 
the softer energy spectra are weaker and shifted to higher frequencies with respect to those 
related to the harder spectra. It leads us to the conclusion that the emission area (CC) does 
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contract when the source undergoes the hard-to-soft state transition. 

Furthermore TLM98 argued that the size of CC strongly depends on the Reynolds (Re) 
number of the accretion flow, and not on just the mass accretion rate. In other words, the 
LHS-HSS transition should be dictated by Re-number values, which is low in the LHS and 
high in the HSS, and the CC size decreases when Re-number increases. TS08 check this 
TLM98 prediction using Cyg X-1 observations by RXTE. First, they find a method to infer 
the Re-number of the accretion flow using the observations. Then they demonstrate that 
when Re-number increases all characteristic frequencies of the power spectrum are indeed 
shifted to higher values as a result of Compton cloud contraction during the LHS-HSS 
transition. 



It is well known from the observations of X-ray BH binaries [see e.g. IWijnands &: van der Klis 



(Il999l ). ST06, ST07 and KVK08] that the power spectra evolve along with the photon spec- 
tra. Here we find that our results for XTE J1650-500 are related to the findings by ST06 for 
Cyg X-1, where the authors, analyzing ~ 10 years of observational data from RXTE archive, 
found strongly correlated characteristics of power spectra and photon spectra (photon index). 
We have to admit that our claims of temporal and spectral correlation in XTE J1650-500 
are based on only the limited set of the BeppoSAX data, while the result s of the extensive 



RXTE data set for XTE J1650-500 bv IShaposhnikov fc Titarchuk I mm . hereafter ST08 



do indicate a strongly correlated temporal and spectral evolution in this source. 

This suggests that the underlying physics of the radiative and oscillatory processes is 
common in X-ray BH binaries. The consistency of the Comptonization model also suggests 
that their spectra emerge from similar Compton cloud geometric configurations whose sizes 
shrink during the state transition. The disk becomes more powerful at softer states and the 
cooling of Compton cloud is mostly dictated by the strong soft disk emission. The photon 
field and surrounding plasma are almost in equilibrium: the plasma temperature is the same 
order of the photon blackbody color temperature (of a few keV). The converging flow located 
in the innermost part of compact cloud is the only place where photons can be upscattered 
to energies of order of mgC^ and higher due to the dynamical Comptonization. 

In the HSS PSs we see the presence of two noise components: VLNF and bknpl (or 
LvLFN and in terms of Lorentzian model). We suggest that the VLFN component is 
associated with the extended disk and bknpl is presumably formed as a result of the diffusive 
propagation of the mass accretion rate perturbation in the compact Compton cloud. In the 
LHS PSs only bknpl is observed, while the VLFN related to the disk component, presumably 
very weak and at very low frequencies, cannot be observed in the BeppoSAX data. 



We also show that in XTE 1650-500 the power of X-ray variability drastically decreases 
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towards the HSS. Note this effect have been known for a number of BHs (see KVK08 and 
TS08 for details of the observations and explanation of this effect, respectively). The relative 
power of the HSS PS with respect of that in the LHS clearly depends on the energy bands. 
We present (for the first time in the literature) energy-dependent PSs of X-ray emission using 
BeppoSAX data. While we find a noticeable dependence of the PS strength on energy, the 
PS shape is found to be almost self-similar in different energy bands (at least for the LHS 
PSs). 

Moreover, we show that the break and QPO low frequencies of the XTE J1650-500 PSs 
which wc derived, correlate with the spectral state, in particular with the PL photon index of 
the energy spectrum. Wc succeed to determine the LF QPO corresponding to two different 
LHS spectra (indices): uqpq = 1.47 ±0.05 Hz, 1.66 ±0.08 Hz related to F = 1.74 ±0.01 and 
1.79 ± 0.06 respectively. Thus our data points based on BeppoSAX data for XTE J1650-500 
are consistent with the index-QPO correlation found by ST08, based on the RXTE data for 
XTE J1650-500. 

We would like to thank anonymous referees for constructive suggestions to significantly 
improve the paper presentation. 
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Table 1: Log of the three observations of XTE J1650-500 



Obs. 


Start time (UT) 


End time (UT) 


LEGS 


MEGS 


HPGSPG 


PDS 








ks 


ks 


ks 


ks 


1 


2001 Sep. 11 10:57:04 


2001 Sep. 12 18:43:06 


20.3 


47.4 


47.0 


22.2 


2 


2001 Sep. 21 18:15:39 


2001 Sep. 23 14:03:40 


15.9 


63.9 


73.3 


30.3 


3 


2001 Oct. 03 16:36:53 


2001 Oct. 04 10:15:56 


7.2 


27.8 


33.0 


12.8 



Table 2. Parameters of the model: WABS x {BBODY + BMC) x SMEDGE x HIGHECUT. HIGHECUT was 
used only for the LHS. Parameters not hsted in the table were held fixed. For HIGHECUT: Ecut — 10 keV; for 
SMEDGE: Eedge — 7.1 keV, index for photo-electric absorption = 2.67 (default value), smearing width = 10 keV. 

Errors are given at the 90% confidence level. 



Model 


Parameters 


1105 


1111 


1215 


1216 




1303 


1308 


WABS 






















n Ro+0.04 




n AC+O.Ol 




0, 


.50 ±0.03 




BBODY 




















kTM (keV) 


u.iy_o.oi 


0.19 ±0.02 


0.36^°:°? 


0.35 ±0.02 


0, 


.34 ± 0.03 


0.33 ±0.02 


BMC 




















r 


1 74+0.01 

^- '^-0.02 


1.81 ±0.02 




2.38 ±0.03 


2, 


.32 ±0.08 


2.34 ±0.04 




kTbb (keV) 


'-'•^'-'-0.01 


0.39 ±0.02 




^■33-0.03 


0, 


r9+0.03 

■J^-0.02 


0.54 ±0.02 




log A 


U.OU_o.o4 


n 20+° °® 

^•^"-"-0.07 


-0 '51+°-i2 


-0 44+°-°^ 




1 06+°-°^ 


-0.831-°^ 




Nsbb/Nbb 


0.8 ±0.4 




1 4+0-9 
■'■•*-0.4 


1 O+0.5 


1, 


1+0.8 
•■'■-0.4 


1 O+0.4 
^^•"5-0.3 


SMEDGE 






















2.04t°;l° 


1 q6+°-i^ 


o r,+0.2 
O.0_o.4 


o e+0.3 
o.0_o.4 


5, 


•)+0.9 


4 r,+0.6 


HIGHECUT 




















Ef (keV) 


114t^ 


114^9 














xVdof 


179.0/164 


167.7/163 


177.6/151 


194.9/153 


164.6/131 


190.7/148 
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Table 3: Parameters of the model: WABS x {BBODY + COMPTT) x SMEDGE. Errors 

are given at the 90% c onfidence level. 

Interval 1105 lUl 

Nn (10^2 cm-2) 0.381^? 0.39^^2 
kT^ (kcV) 0.30 ± 0.01 0.31 ± 0.01 
TO (keV) 0.45 ± 0.03 0.46i[J:|]^ 

k% (keV) 311? 321^ 

r 2.15 ±0.03 2.12 ±0.06 

SMEDGE-Tmax 1-741^2 1-7 ± 0.2 

xVdof 237.0/165 210.4/164 



Table 4. Characteristics of BH XTE J1650-500 variability components of the power spectrum for Low Hard State 
(TOO-1) and High Soft state (TOO-2) [see Fig. E] (bottom panel) for the related power spectra]. The model is a 
superposition of two broad Lorentzians (L^, L/^) plus one or two narrow Lorentzians (Lgpoi, Lqpo2)- Parameters in 
square parentheses are frozen in the fits. Errors are given at the 90% confidence level. See text about the 

identification of the component. 



Components 0.1-2 keV 2-10 keV 



Hard State (TOOl) 






Lb, l^b max (Hz) 


0.591°:?^ 


n Qo+0.08 


L/i, t'/i max (Hz)" 


23 ±2 


2 

^■^-0.3 


VLNF, VVLNF max (Hz) 






Lqpoi, ^qpo (Hz) 


1.55 ±0.08 


1.53 ±0.04 


LqP02, t'QPO (Hz) 




3.08 ±0.07 


Soft State (T002) 






Lfo, l^b max (Hz) 


36 ±3 


[35] 


L/i, t'/i max (Hz) 






VLNF, VVLNF max (Hz) 


(5^^) X 10-2 


{A.Qt\l) X 10-^ 



'^See section 4.1, for details about the identification of the 
Lh component in the 0.1-2 and 2-10 keV energy intervals 



5. Energy dependence of PS component characteristics in Hard state (TOO-1) for the Lorentzian model. 
Parameters in square parentheses are frozen in the fits. Errors are given at the 90% confidence level. 



Components 


0.1-2 keV 


2-4 keV 


4-10 keV 


8-15 keV 


15-200 keV 


Lb, l'b,m.ax (Hz) 
Lh, t^h,max (Hz)" 

Lqpoi, i^q (Hz) 
Lgpo2, 1^0 (Hz) 


0.59tn^ 
23 ±2 


0.38 ± 0.06 
2.1 ±0.2 


n QR+0.08 
0.36Io.o4 

^■3-0.5 
o no+O.lO 

O.UO_Q gg 


0.38tn^ 

o tr+1.4 

[1.54] 
[3.03] 


q q+1.0 
> 50 



"^See section 4.1 for details about the identification of the Lh component in the 
0.1-2, 2-4, 4-10, 8-15 and 15-200 keV intervals 



Table 6. Best model fit parameters used to fit the power spectra of TOO-1 in different energy bands. The model: 
BKNPL + LORENTZlANl + LORENTZIAN2. Parameters in square parentheses are frozen in the fits. Errors are given at 

the 90% confidence level. 



Parameter 


0.1-2 keV 


2-4 keV 


4-10 keV 


2-10 keV 


8-15 keV 


15-200 keV 


ai 




[0.] 


[0.] 


[0.] 


[0.] 


[0.] 


[0.] 






0.7 ±0.1 


1.20 ±0.02 


1.05 ±0.04 


1.098 ±0.013 


0.96 ±0.03 


0.49 ±0.10 


l/br (Hz) 




0.28 ±0.03 


0.345 ±0.016 


0.36 ±0.03 


0.35 ±0.01 


0.37 ±0.04 


0.4 ±0.4 


uli (Hz) 






1.49 ±0.09 


1.52 ±0.04 


1.53 ±0.03 


1.54 ±0.04 


1.54 ±0.15 


FWHMli 


(Hz) 




1.1 ±0.3 


0.8 ±0.3 


0.9 ±0.2 


0.46 ±0.18 


1.1 ± 1.4 


i^Li (Hz) 








[2j^li] 






FWHMl2 


(Hz) 




1.6 ±0.8 


0.9 ±0.3 


0.95 ±0.28 


1.1 ±0.5 


1±1 


xVdof 




76.6/86 


144.1/140 


117.8/106 


192.6/182 


95.9/68 


8.6/13 
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Table 7. Best fit parameters of the models used to fit the power spectra in two different 

time intervals (I101-I107, I109-I112) of TOO-1 in 4-10 keV . Model: bknpl + 
LORENTZlANl + LORENTZIAN2. Parameters in square parentheses are frozen in the fits. 

Errors are given at the 90% confidence level. 



Parameter 


I101-I107 


I109-I112 


CCi 




[0.] 


[0.] 


Ot2 




0.961 ±0.025 


1.23 ±0.12 


Uhr (Hz) 




0.35 ±0.03 


0.37 ±0.07 


i^Li (Hz) 




1.47 ±0.03 


1.66 ±0.08 


FWHMii 


(Hz) 


0.63 ± 0.25 


1.2 ±0.7 


J^L2 (Hz) 






[2^i] 


FWHMi2 


(Hz) 


0.5 ±0.2 


0.9 ±0.6 


xVdof 




110.4/107 


32.6/42 



Table 8. Best fit parameters of the models used to describe the power spectra of TOO-2 (HSS). In 0.1-2 keV, two 
models: 1st model: PL + bknpl; 2nd model: PL + constant. The first index of the bknpl was fixed to 0. In 2-10 

keV. Errors are given at the 90% confidence level. 



Model 



Parameter 



0.1-2 keV 



2-10 keV 



PL+BKNPL 



PL+CONST 



a 

NPL 
«! 

0.2 

Ubr (Hz) 

Nbknpl 
a 

NpL 

CONST 
xVdof 



1.0 ±0.3 
(4.011°) X 10-4 
[0.] 



(2.0l°i) X 10-' 



126/102 



2 ± 1 

(3^r) X 10-8 

(1.29 ±0.35) X 10-5 
16.6/20 
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Fig. 1. — Light curves of XTE J1650-500 in two energy bands: 2-3 keV {top), 20-60 keV 
{bottom). 
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Fig. 2. — Ratio between MECS mean count rate in the time interval 1202 (TOO-2) and that 
in interval 1207 (TOO-2) as a function of energy. 




' JH. ■ ■ I 1 1 I I 1 1 1 I 
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Fig. 3. — Evolution of the EF{E) energy spectra with time. The LHS spectra (1105 and 
nil) show the key role of the thermal Comptonization of soft (disk) photons in the extended 
hot Compton cloud (corona), while the HSS spectra (I215-I216, 1303, and 1308) show that 
prominent contribution of the soft blackbody-like component due to the direct emission 
from the disk. In this case the hard X-ray component (i.e., the steep PL) is formed in the 
relatively cold (temperature of a few keV) compact region which is presumably a converging 
flow into the BH. 



- 35 - 
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0.1 



' r 

\ 
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0.15 



0.2 



BMC normalization l^g/d 



Fig. 4. — Photon index evolution as a function of BMC normalization L^g/df^^ (which is 
proportional to mass accretion rate in the disk). Colors correspond to TOOl(black), T002 
(red), T003 (green). 
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1 1 1 1 — I — I — r-r-| 1 1 1 1 — I — I — r-r-| 1 1 1 1 — i — i — r-r-| : 

1 1 1 — I I I I I I 1 1 1 — I I I I I I 1 1 1 — I I I I I I E 

— I — I — I — I 1 1 1 1 1 1 — I — I — I 1 1 1 1 1 1 — I — I — I 1 1 1 1 1 
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1 1 1 — I I I I I I 1 1 1 — I I I I I I 1 1 1 — I I I I I I E 

— I — I — I — 1 1 1 1 1 1 — I — I — I — 1 1 1 1 1 1 — I — I — I — 1 1 1 1 1 1 — \ 

+++ V ^ 

I I I I I I I I I I I I I I I I I I I I I I I I I I I IE 

0.1 1 10 100 

channel energy (keV) 

Fig. 5. — Residuals (in standard deviation units) of the source spectrum in the time interval 
1108 using different best-fit spectral models in different energy bands (see text). From top 
to bottom. Panel 1: XSPEC model WABS x {BB + BMC) x HIGHENCUT fits to data 
in the 0.1-200 keV range, xV(dof)=430/165. Panel 2: our best-fit model {WABS x {BB + 
BMC) X SMEDCE x HIGHENCUT) fits to the same data, xV(dof)=150/164. Panel 3: 
XSPEC model WABS x{BB + BMC + LAOR)x HIGHENCUT, xV(dof)=300/163 fits to 
the data in the 0.1-200 keV. Panel 4'- The same spectral model, adopted to fit the spectrum 
of panel 1, is used to fit the spectrum in the 1.5-60 keV, x^/(dof)=200/lll. Panel 5: The 
same spectral model, adopted to fit the spectrum of panel 3, is used to fit the spectrum in 
the 1.5-60 keV energy range, x^/(dof) = 136/109. Panel 6: The same model and parameters 
of panel 5, is used to fit the 0.1-200 keV energy range, x^/(dof) = 1315/164. 
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10 ^ 0.01 0.1 1 10 0.01 0.1 1 10 100 



Frequency (Hz) Frequency (Hz) 

Fig. 6. — The LHS PS (blue) vs. HSS PS (red) for two different energy bands: 2-10 keV 
{left panels) and 0.1-2 keV {right panels). In tfie 2-10 keV, top panels: the HSS PS is fitted 
by the sum of a PL plus a constant (see Table [8]), while the LHS PS is fitted by a bknpl 
plus two Lorentzians (see Fig. 5); bottom panel, frequency x power diagram: the HSS PS is 
fitted by the sum of two broad Lorentzians, while the LHS PS is fitted by the sum of 2 broad 
Lorentzians plus two narrow Lorentzians. It is apparent that the relative power of the HSS 
PS (with respect to the LHS PS) strongly depends on the energy band. 
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Fig. 7. — Top panel: the LHS power spectra in different energy bands (see bottom panel for 
the correspondence between color and energy intervals). The decrease of the variability power 
with energy is evident from this Figure. Bottom panel: the LHS power spectra as above, 
shifted, for clarity of display, along the ordinate direction by different factors, depending on 
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Frequency (Hz) 



Fig. 8. — The LHS power spectra of XTE J1650-500 in the time intervals from 1101 to 1107 
and from 1109 to 1112 (see text for details). The power spectrum (red line) corresponding 
to the softer energy spectrum is shifted to higher frequencies with respect to that (blue line) 
corresponding to the hard energy spectrum (see related energy spectra in Fig. [3]). The fit 
with a BKNPL plus two Lorentzians is also shown. Best fit details in Table [61 
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Fig. 9. — Evolution of the Compton cloud area during the spectral transition exemplified 
by the progressive number of the time intervals from TOO-1 to TOO-3 (see text). Colors 
correspond to TOOl (black), T002 (red), T003 (green). 
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Fig. 10. — Evolution of the Comptonization fraction / = A/{l + A) during the spectral tran- 
sition from LHS to HSS. Colors correspond to TOOl (black), T002 (red), T003 (green). 



